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Comment on "Generalized Dynamic Scaling for Critical Relaxation" 
[Phys. Rev. Lett. 77, 679 (1996)] 

In the work by Zheng M the temporal evolution of a spin system that is 
quenched from a high-temperature initial state to the critical point T c is 
studied by means of scaling analyses and Monte Carlo simulations. Special 
attention is paid to the time dependence and scaling behavior of the order- 
parameter m(t), when the system - the concrete case studied was the 3-state 
. Potts model - is quenched from a state with a finite initial magnetization 

\ The analysis of Ref. is built upon previous work [g, |J where the 

' influence of a small mo on relaxation processes in finite systems was studied. 

When, after the quench, the time-dependent correlation length ~ t l l z 
(where z denotes the dynamic equilibrium exponent) has grown much larger 
CN '■ than macroscopic scales, e.g. the lattice spacing a, the order parameter 

m(t, m ) should, due to the loss of intrinsic scales, exhibit universal behavior. 
Especially for mo — > in a system of finite linear size L (3> a) the scale 
invariance can be expressed in the form 

m(t, L, m ) = b- p t" ' mitb-^Lb- 1 , m b xo ) . (1) 

S ■ 

The exponent Xq in this relation, the scaling dimension of the initial mag- 
netization, is an independent critical exponent that in general is different 
from the scaling dimension /3/u of the equilibrium magnetization || |2|. By 
mapping Monte Carlo profiles from systems of different size onto each other 
. ^ ' and by employing Eq. (P , it was possible to estimate the new exponent Xq 

X': II- 

For larger m , Eq. fll]) as it stands becomes certainly wrong since mo 
saturates. Nevertheless one would expect that in a stage where £(t) has 
grown much larger than a the temporal evolution of m(t) is still universal 
and should exhibit scaling. To verify this, Zheng allows in the last argument 
on the right-hand side of ([!]) for a general function m (mo). Similar to the 
notation for small mo, this function is written in the form m = \fo( b > m o) mo, 
and xo(b, m ) is determined as a function of m . The results presented in 
Figs. 1,2, and 3 of Ref. prove that scaling occurs for general m , and there 
is no doubt on the correctness of these results. 

What we want to point out in this Comment is the following: In order 
to take into account the initial state in the scaling analysis the field that has 
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to enter the description is actually the magnetic field Hq that is imposed in 
the initial state (and switched off at t — 0) and not the initial magnetiza- 
tion. Only for small H one can effectively replace H by mo for the initial 
state is paramagnetic and responds linearly on small magnetic fields. As a 
consequence, it was not necessary to distinguish between H and m in most 
of the earlier work|FJ, |3|. Now, the original Xq in (|l|) is actually the scaling 
dimension of H Q - we call it i/q in the following -, and the distinction becomes 
important in a regime where the response is not linear anymore. As a result, 
Eq. 1 should be replaced by 

m{t, L, H ) = b- p/v m(t h~ z , L &~\ H b yo ) , (2) 

and especially yo should be simply a constant for all Hq. Assuming the latter 
and using that itlq = tanh(ifo) m our high-temperature initial state, the 
function xo(b, m ) occurring in ([[]) can even be calculated. The result is 
given by 

Xo(b, m ) = - — - {log [tanh (b yo artanhm )] — logm } . (3) 

To corroborate this by Monte Carlo simulations, we calculated a number 
of magnetization profiles for the d = 2 Ising model for lattices with sizes 40 2 
and 80 2 , such that 6 = 2, and with periodic boundaries. Further we used the 
heat-bath algorithm and sequential updating. In order to rescale the profiles 
we adopted the literature value z = 2.170 and used the exact exponent 
(3/v = 1/8. In Fig. 1 we show yo (squares) and xq (full circles) for b = 2. 
The errors were estimated. The value Xq = y$ ~ 0.54 for m = is derived 
from the short-time exponent 9 = (xq — (3/v)/ z ~ 0.19 [f|, |(|. Further, we 
have obviously Xq = for mo = 1, whereas yo can not be determined at this 
point since H = oo. For comparison we also show the result of Eq. (|3|) for x 
(solid line), which is consistent with the numerical data. Most importantly, 
our results for y are consistent with a constant (dashed line) for all H . 

Eq. (^|) can straightforwardly be generalized for the p-state Potts model. 
The result for p = 3 with xq(2, 0) = 0.3 is also shown in Fig. 1 (dotted line). 
The curve is qualitatively consistent with the data of Ref . , however, it does 
not provide a quantitatively correct fit. This discrepancy could be related to 
the fact that in Zheng's simulations the initial state is not a thermal equilib- 
rium state but has a sharp m . We do not believe that for the general Potts 
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Figure 1: The scaling dimensions xo (full circles) and yo (squares) as a func- 
tion of mo from lattices of size 40 2 and 80 2 . The solid line represents the 
result of Eq. (3) for xq(2, mo), and the dotted line is the corresponding curve 
for the 3-state Potts model. 

model (|2|) fails. Our simulations always started from a thermal equilibrium 

state. 

Let us finally remark that a situation in many respects analogous to the 
one discussed above is met in equilibrium systems with a planar surface (semi- 
infinite geometry) at bulk criticality. To the high-temperature initial state 
in the dynamical problem corresponds a surface that strongly suppresses the 
order, and the analogue to H is a magnetic field Hi that is restricted to a 
microscopic region near the surface. As discussed in Ref. ^ it is of particular 
importance to keep fields and densities distinguished in the context of surface 
critical phenomena, and much can be learned from there for nonequilibrium 
critical dynamics (and vice versa). 
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